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Abstract 

Various lepton-flavor violating (LFV) processes in the supersymmetric standard 
model with right-handed neutrino supermultiplets are investigated in detail. It is 
shown that large LFV rates are obtained when tan (3 is large. In the case where 
the mixing matrix in the lepton sector has a similar structure as the Kobayashi- 
Maskawa matrix and the third-generation Yukawa coupling is as large as that of 
the top quark, the branching ratios can be as large as Br{fi — > 67) ~ 10~^^ and 
Br[T jjL^) ~ 10^'', which are within the reach of future experiments. If we 
assume a large mixing angle solution to the atmospheric neutrino problem, rate for 
the process r fij becomes larger. We also discuss the difference between our case 
and the case of the minimal SU (5) grand unified theory. 



1 Introduction 



Lepton-flavor violation (LFV), if observed in a future experiment, is an evidence of new 
physics beyond the standard model, because the lepton-flavor number is conserved in 
the standard model. Since the processes do not suffer from a large ambiguity due to 
the hadronic matrix elements, detailed analysis of the LFV processes will reveal some 
properties of the high-energy physics. 

One of the minimal extensions of the standard model with LFV is the model with 
non- vanishing neutrino masses. If the masses of the neutrinos are induced by the seesaw 
mechanism [1] , one has a new set of Yukawa couplings involving the right-handed neutri- 
nos. Introduction of the new Yukawa couplings generally gives rise to the flavor violation 
in the lepton sector, similar to its quark sector counterparts. In non-supersymmetric stan- 
dard models, however, the amplitudes of the LFV processes are proportional to inverse 
powers of the right-handed neutrino mass scale which is typically much higher than the 
electroweak scale, and as a consequence such rates are highly suppressed. 

If the model is supersymmetrized, the situation becomes quite different. LFV in 
the right-handed neutrino Yukawa couplings leads to LFV in slepton masses through 
renormalization-group effects [2] . Then the LFV processes are only suppressed by powers 
of supersymmetry (SUSY) breaking scale which is assumed to be at the electroweak scale. 
Especially, in a previous paper [3], we pointed out that a large left-right mixing of the 
slepton masses greatly enhances the rates for the LFV processes such as jj ^ and 
r fi'-f. Due to this effect, they can be within the reach of near future experiments even 
if the mixing angle of the lepton sector is as small as that of the quark sector. 

In this paper, we will extend the previous analysis. We are interested in the following 
processes, 

• r ^ i^-f, 

• II ^ eee, 

• ii-e conversion in nuclei, 

and calculate formulas for the interaction rates of the above processes. In our calculation. 
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we fully incorporate the mixing of the slepton masses as well as the mixings in the neu- 
tralino and chargino sectors. Also the lepton Yukawa couplings in higgsino-lepton-slepton 
vertices are retained, which yield another type of enhanced diagrams in the large tan (5 
region. Then we will discuss how large the interaction rates can be, assuming the radia- 
tive electroweak symmetry breaking scenario [4] . We find that a large value of tan (3 is 
realized with relatively light supcrparticle mass spectrum, and thus the interaction rates 
can indeed be enhanced. For the right-handed neutrino sector, wc will mainly consider 
the case where the Yukawa couplings of the right-handed neutrinos arc similar to those 
of the up-type quarks. We will also discuss the case of large mixing between the sec- 
ond and third generations, suggested by atmospheric neutrino problem. In our numerical 
analysis, we impose the constraints from the negative searches for the SUSY particles, as 
well as the constraint from the muon anomalous-magnetic dipole-moment g — 2 to which 
supcrparticle loops give non-negligible contributions especially in the large tan (5 region. 

The organization of our paper is as follows. In the subsequent section, we will review 
LFV in slepton masses in the presence of the right-handed neutrinos. In Section 3, we 
will give formulas of the interaction rates of the various LFV processes. Results of our 
numerical study are given in Section 4. In Section 5, after summarizing our results, we will 
compare our case with the case of the SU (5) grand unification briefly. Renormalization- 
group equations relevant to our analysis are summarized in Appendix A. In Appendix B, 
we describe the interactions among neutralino (chargino)-fermion-sfermion. In Appendix 
C, we will give formulas of the SUSY contribution to g — 2. 

2 LFV in scalar lepton masses 

Throughout this paper, we consider the minimal SUSY standard model (MSSM) plus three 
generation right-handed neutrinos. In this case, the superpotential is given by 

+fxe^pH?H^ + \M'jN^N'r, (1) 

where Lj represents the chiral multiplet of a SU (2) l doublet lepton, E^ a SU (2) ^ singlet 
charged lepton, Nf a right-handed neutrino which is singlet under the standard-model 
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gauge group, Hi and H2 two Higgs doublets with opposite hypercharge. Similarly Q, U 
and D represent chiral multiplets of quarks of a SU{2)l doublet and two singlets with 
different U{1)y charges. Three generations of leptons and quarks are assumed and thus 
the subscripts i and j run over 1 to 3. The symbol is an anti-symmetric tensor with 
ei2 — 1- The Yukawa interactions are derived from the superpotential via 



1 d'^W , , ^ 



(2) 



SUSY is softly broken in our model. The general soft SUSY breaking terms are given 



as 



l^'^Ml^u^^ + h.c.) 
+ImiBIbI + -M2WIWI + -M^G^G" + h.c). 



(3) 



2 2 2 

Here the first four lines are soft terms for sleptons, squarks and the Higgs bosons, while 
the last line gives gaugino mass terms. 

We now discuss LFV in the Yukawa couphngs. Suppose that the Yukawa coupling 
matrix //-' and the mass matrix of the right-handed neutrinos Ml^ are diagonalized as 
fiiS'^^ and MjiiS'^^ , respectively.^ Then, in this basis, the neutrino Yukawa couplings pj 
are not generally diagonal, giving rise to LFV. An immediate consequence is neutrino 
oscillation. Writing fl^ = U'^ fvkV^^ with [/, V unitary matrices, we obtain the neutrino 
mass matrix induced by the seesaw mechanism 

.,2 



fjM-'U X y sin' P 



\ 



( 



1/2 



1/3 



Mr2 



f: 



ul 



u 



1/2 



Mrs J 



Us J 



V 



X — sm 3, 
2 



(4) 



^We can always choose and Ml^ to be diagonal by using unitary transformations of L, and 
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where ^v^ = {hi)^ + (/ia)^ ^ (174GeV)2 and tan/? = {h2)/{hi). (Here, (• • •) stands for 



the vacuum expectation value of the quantity.) Throughout this paper, we assume that 
Ml, is proportional to the unit matrix Ml^ — MrS'^^ , for simphcity. Then, if we disregard 
possible complex phases in [/, the above can be rewritten as 



Thus as far as ^ 7^ 1 and the mass eigenvalues are non-degenerate, we have neutrino 
oscillation which is a target of current and future experiments. 

The smallness of the neutrino masses implies that the scale Mr is very high, ~ 10^^ 
GeV or even higher. In the standard model with right-handed neutrinos, the flavor 
violating processes such as — > 67, t ^ /I'j etc., whose rates are proportional to inverse 
powers of M^j, would be highly suppressed with such a large Mji scale, and hence those 
would never be seen experimentally. 

However, if there exists SUSY broken at the electroweak scale, we may expect that 
the rates of these LFV processes will be much larger than the non-supersymmetric case. 
The point is that the lepton-fiavor conservation is not a consequence of the standard- 
model gauge symmetry and renormalizability in the supersymmetric case, even in the 
absence of the right-handed neutrinos. Indeed, slepton mass terms can violate the lepton- 
fiavor conservation in a manner consistent with the gauge symmetry. Thus the scale 
of LFV can be identified with the electroweak scale, much lower than the right-handed 
neutrino scale Mr. However, an order-of-unity violation of the lepton-fiavor conservation 
at the electroweak scale would cause disastrously large rates for // — > 67 and others. 
Also, arbitrary squark masses result in too large rates for various flavor-changing- neutral- 
current processes involving squark loops. To avoid these problems, one often considers 
that the sleptons and the squarks are degenerate in masses among those with the same 
gauge quantum numbers in the tree-level Lagrangian at a certain renormalization scale. 
In the following, we will assume a somewhat stronger hypothesis that all SUSY breaking 
scalar masses are universal at the gravitational scale M = mpi/y/Sn ~ 2 x lO^^GeV, i.e., 
we adopt the minimal supergravity type boundary conditions. Thus we will consider the 
following type of soft terms. 




V 9 

V X — sin 
2 



(5) 
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• universal scalar mass (mo), 

all scalar masses of the type (mj)] and ml. {i — 1,2) take common value ml, 

• universal 74-parameter, 

A^j — afJmQ with a being a constant of order unity, 

at the renormalization scale M.^ As for the gaugino masses, for simplicity, we choose 
the boundary condition so that they satisfy the so-called grand unified theory (GUT) 
relation at low energies. Note that the universal scalar masses are given in a certain 
class of supergravity models with hidden sector SUSY breaking [5]. Those soft SUSY 
breaking terms suffer from renormalization via gauge and Yukawa interactions, which can 
be conveniently expressed in terms of the renormalization-group equations (RGEs). The 
RGEs relevant in our analysis will be given in Appendix A. An important point is that, 
through this renormalization effect, LFV in the Yukawa couplings induces LFV in the 
slepton masses at low energies even if the scalar masses are universal at high energy. Due 
to this fact, lepton- flavor conservation is violated at low energies. 

We can solve the RGEs numerically with the boundary conditions given above. It 
is, however, instructive to consider here a simple approximation to estimate the LFV 
contribution to the slepton masses. Since the SU{2)l doublet lepton multiplets have the 
lepton-flavor violating Yukawa couplings with the right-handed neutrino multiplets, the 
LFV effect most directly appears in the mass matrix of the doublet sleptons. The RGEs 
for them can be written as (see Appendix A) 



1 

+ 



mi/^/. + /^/.m|)i + 2(/>?/. + m^^/t/. + 4^.)] . (6) 



167r2 

Here (/i^(m|)])]v[ssM denotes the RGE in case of the MSSM, and the terms explicitly 
written are additional contributions by the right-handed neutrino Yukawa couplings. An 
iteration gives an approximate solution for the additional contributions to the mass terms 

(Ami); « -^^^ ^^ {SmUmi + 2(AlA.)i) 

^In fact, there is another SUSY breaking parameter B, which gives a mixing term of the two Higgs 
bosons hi and /i2. For a given value of tan/3, we fix this parameter B (and also the SUSY invariant 
Higgs mass fi) so that the Higgs bosons have correct vacuum expectation values, {hi) = vcos(}/\/2 and 
(/12) = vsml3/s/2. 
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- _!!^(6^2a>M/„)|. (7) 
where we have used the universal scalar mass and A-parameter conditions. In Eq. (7), 

{flf^i = fLft' = v,:\Wv'^, (8) 

so that the slepton mass indeed has the generation mixing if V differs from the 

unit matrix in the basis that the charged lepton Yukawa coupling fi are diagonal. 

Lack of our knowledge on the neutrino Yukawa couplings prevents us from giving a 
definite prediction of the slepton mass matrix, and thus the rates of the LFV processes. 
Nevertheless, it is important to study how large the interaction rates for the LFV pro- 
cesses can be for some typical cases and to see whether those signals can be tested by 
experiments. In this paper, we shall consider the following typical two cases: case 1 ) the 
mixing matrix V is identical to the Kobayashi-Maskawa (KM) matrix in the quark sector 
Vxivh case 2) the mixing matrix is given so that it can explain atmospheric neutrino 
deficit by the large-mixing v^-v^ oscillation. In the latter case, we only consider r — > /X7, 
the generation mixing between the second and third ones. 

3 Interaction rates for LFV processes 

In this section we give formulas of the interaction rates for the LFV processes we consider. 
Results of our numerical calculation will be given in the next section. 

We first explain how the rates for — > 67 and r — > //7 can be enhanced compared 
with the naive expectation when tan (5 is large. Here, we consider in the basis where the 
neutralino/chargino interactions to the leptons and the sleptons are fiavor diagonal and 
the effect of the fiavor violation in the lepton sector is involved by the mass insertions 
(m|)i {i 7^ j). First, let us consider contribution from winos and bino, the SU {2)lxU (l)y 
gauginos, neglecting the mixing in the chargino/neutralino sector. A naive estimate on 
the branching ratio yields 

5r(/,-/.7)oc^^^-^, (9) 

where 777.5 is the typical mass of superparticles, a the fine structure constant and Gp the 
Fermi constant. The contribution from Feynman diagrams Fig. 1(a) and (b) follows this 
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estimate. However, as emphasized in our previous paper [3], the diagram Fig. 1(c) which 
picks up the left-right mixing of the sleptons and exchanges the bino in the loop can 
give much larger contribution, when /i tan /3 is much larger than the masses of the other 
superparticles. Indeed we estimate the ratio of the amplitudes 



Amp. ( 1 . c) Mimlj^jj tan /? Mi 



(10) 



Amp.{l.a + b) mi.mg ms ms^ 

with mi- being the charged lepton Ij mass. In Ref. [3], we numerically showed that this 
enhancement really occurs for the case of large /i tan /3. 

If we take account of the gaugino-higgsino mixing in the chargino/neutralino sector, 
we find another type of diagram which enhances the amplitude when tan /3 is large but 
IJ, is comparable to the masses of the other superparticles. It is shown in Fig. 2. In this 
diagram, one has the mixing between the higgsino and the gaugino which is proportional to 
V sin (3, the vacuum expectation value of h2, and involves the Yukawa coupling of higgsino- 
lepton-slepton, fi. = —\/2mi./{vcosP). The sleptons inside a loop are left-handed ones. 
Thus the amplitude is proportional to tan and 

Amp.{2) 

—r^tan^. (11) 

Note that this type of diagram includes neutralino-exchange graphs as well as a chargino- 
exchange graph. 

In this work, we are interested in the following LFV processes; — 67 and r /i7, 
ix~ — > e~e~e^ and ^i-e conversion in nuclei. To obtain the interaction rates for these 
processes, we perform full diagonalization of the slepton mass matrices numerically and 
consider mixing in the chargino and neutralino sectors. 

We write the interaction Lagrangian of fermion-sfermion-neutralino as 

>C = K^TxPr + nISPl)x\]x + h.c. (12) 

In this section, /j (/ = IjUjdjU) represents a fermion in mass eigenstate with the gener- 
ation index i {i = 1,2,3), and fx a sfermion in mass eigenstate. The subscript X runs 
from 1 to 3 for u and from 1 to 6 for the other sfermions, /, d and u. A neutralino mass 
eigenstate is denoted by {A — 1, • • • , 4) and Pr^l = |(1 i 75). The coefficients 
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and nI^x depend on the mixing matrices of the neutrahno sector and of the sfermions. 
Their exphcit forms will be given in Appendix B. Similarly the fermion-sfermion-chargino 
interaction is written as 

+u,{Cfi'iPR + Cl:i'iPL)x\dx + h.c, (13) 

where (A = 1, 2) is a chargino mass eigenstate. The explicit forms of the coefficients 
can also be found in the Appendix B. 

3.1 Effective Lagrangian for LFV processes 

As a first step to compute the LFV rates, let us write down the effective interactions (or 
amplitudes) relevant for our purpose. 

3.1.1 IJ^l-i* 

The off-shell amphtude for Ij — > Z~7* is generally written as 

T = ee"*M,(p - q) [g^a(Af + A^Pr) + m,^.za„^/(yl^Pi + A^Pr)] u,{p), (14) 

in the limit of g — > with q being the photon momentum. Here, e is the electric charge, e* 
the photon polarization vector, Ui (and Vi in the expressions below) the wave function for 
(anti-) lepton, and p the momentum of the particle Ij. In the present case, the Feynman 
diagrams contributing to the above amplitude are depicted by Fig. 3. Each coefficients in 
the above can be written as a sum of the two terms, 

where A^^^'^ and A^^^'^ stand for the contributions from the neutralino loops and from 
the chargino loops, respectively. We calculate them and find that the neutralino contri- 
butions are given by 

.(n)L _ 1 ^,R(0^,R(0' 1 1 

^1 — r^c_9^^iAX^^l 



5767r2 ^'^^ m? (1 - xaxY 

Ix 
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,{n)L 



x(2 - 9xAx + iSx\x - IIJ'Iy 

1 1 r L(Z) L(/)* 1 



327r2 



m 



6(1 - xyix^ 



X (1 - 6xax + ^x\x + 2a;^x - ^^ax ^ ^ax) 



a 



1,2), 



(15) 



(16) 
(17) 



where xax = /ni^ is the ratio of the neutrahno mass squared, M?o , to the charged 



slepton mass squared, ^f^- (Summation over the indices A and X are assumed to be 
understood.) The chargino contributions are 



a(c)l 

^2 



1 ^R{l)^R{l)* 1 



5767,2 -^A^-. AX ^2^(l_^^^)4 
X {l6 - 45a; AX + 36xax - ^x^x + 6(2 - Sxax) Inx^x} , 



(18) 



1 



327r2 m? 



^^^^^^■^^ 6(1 - XAxY 



X (2 + 3xax - ^x\x + 3:^Ax + 6a:; ax In xax) 
M— 1 

+^fiS^^S*^7:^ ^(-3 + ^XAx - x\x - 21nxAx) 

mi. (1 - xaxY 

A^^%^R (a = 1,2). 



(19) 
(20) 



Here, xax — M?- / , where M~- and m^^ are the masses for the chargino the 
sneutrino z/x, respectively. 



3.1.2 IJ 



Ia Ia Ia 



We next consider the process l~ — > If (including e~e~e+). The effective 

amplitude consists of the contributions from the Penguin-type diagrams and from the 
box- type diagrams. The former contribution can be computed using Eq. (14), with the 
result 



T 

'y— penguin 



Uiipl) q'lMiPL + A^Pr) + mi.ia^pq^{A^PL + A^Pr) Ujip) 



\R, 



X^Ui{p2)'y"Vi{ps) - (pi^P2)- 



(21) 
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Furthermore, there are the other Penguin-type diagrams in which the Z boson is exchanged 
as shown in Fig. 4. This amphtude is 

Tz-penguin = ^Ui{p^)r {FlPl + FRPR)uj{p)Ui{p2)riZ[PL + Z'^PR)Vi{p^) 

rriz 

-(pi^P2), (22) 

where F^r) — F^j^-j + F^^^ . The chargino contribution F^^^ and the neutrahno contri- 
bution are^ 



^ 167r2 



R(l)^R{l)* r, 

1/T-il^nll/l-ll T-tn II t T I An \ I It I 1~m 

(23) 



^{x,A,B) 7. '^{XA,B) 



Ft' = 0, (24) 

i-,(n) ^iAX^jBX* {On)a3{On)b3 " (^Af) A4(QAr)i?4 ^ \ .^ s 

= n {J'{X,A,B)+^<^(X,A,B) , (25) 



167r2 

^(n) ^ (26) 

Here, 0^, ^ and Ojv are orthogonal matrices to diagonahze the mass matrices of the 
chargino and neutrahno (see Appendix B), and F(^x,a,b) and G(x,a,b) are given by 

F(x,A,B) = mxAxH [— z , (27) 

xax - xbx V 1 ~ ^Ax 1 - a^^BX y 

(^{x,A,B) = —-2 — -j ■ (28) 

TTT'i^ Xax — Xbx \ ^ — xax ^ — xsx J 

In these functions, M^^ and mf^ denote neutrahno mass and charged slepton mass in the 
neutrahno contribution, and chargino mass and sneutrino mass in the chargino contribu- 
tion. And, in Eq. (22) the coefficient Z^^^^ denotes Z boson couphng to charged lepton 
Il(R), that is, 

^i(R)=r3L(H)-QLsm'^iy, (29) 

where T^l^r) and represent weak isospin (Tg^ — ~h ~ ^) electric charge 
(Qim = -1) of Il{r) respectively. 

^The Penguin-type diagrams of Z boson contributing to the LFV events do not necessarily need to 
have chirahty flip of lepton as /i — > ej. Therefore, the diagrams picking up Yukawa coupling of higgsino- 
fermion-sfermion can not become the dominant contribution in Z boson Penguin-type diagrams and we 
neglect them in the above equations. 
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The box-type Feynman diagrams are given in Fig. 5, and we can write their amphtude 

as 

Tbox = B^e"^ Ui(j)i)'y"'PLUj{p) Ui{p2)laPLVi{p3) 

+Sfe^ Ui{pi)-f'^PRUj{p) Ui{p2)-faPRVi{p3) 
+S^e^ {Ui{pi)-f°'PLUj{p) Ui{p2)-faPRVi{p3) " (Pl^P2)} 
+5^e^ {Ui{pi)-f°'PBUj{p) Ui{p2haPLVi{p3) - {pi^P2)} 
+B^e^ {Ui{pi)PLUj{p) Ui{p2)PLVi{p3) - (Pl^P2)} 
+3^6^ {Ui{pi)PRUj{p) Ui{p2)PRVi{p3) - (pi^P2)} 
+^46^ {Ui{pi)cr^i,PLUj{p) Ui{p2)(j'"'PLVi{p3) - {pi^P2)} 

+Bfe^ {ui{pi)a^,PRUj{p) Ui{p2)a^'' PRVi{p3) - {pi^P2)} , (30) 



where 



^L,R ^ ^{n)L,R ^ ^(c)L,R ^ 1^ . . . ^ 4). (31) 



The first term represents the neutrahno contribution, which we find to be 

2 rj(n)i _ 1 T isrR{l)* atR{1) ArR(0* /\r-R(0 



+/4(w,y)M,.,M,o^iV/iJ*iv5«*iv5^ (32) 

^ T . J A7^('> A7^(') A7^^'> A7^(') ^ A/^^')* M^^^> A7^(') A7^(') 

^^4(A,B,X,y) ^^iAY^^iBY iBX + -'^'AX ^^jAF ^^iBY^^iBX 

'^^jAX ^^lAY ^^iBY^^iBX] 



\hiA,B,x,Y)M^^^M^o^Nf2i*N^3N^ii'N^^^^^ (33) 
2 rj(n)L _ . A/f „ A/f . / /xr-^C)* /vr-^^') /\r^(0* Ar^O) 



+^<i?ArfJ*Arf«Arfi$}, (34) 

e^sf^^ = ^/4(AB,x,y)M,o^M,o^iVfiJ*ArS*ArS?Arfj^ (35) 
^(n)i? ^ sW^I^^^ (a=l,---,4). (36) 



The chargino contribution is 



^2o(c)i:' _ 1 T ^R{l)*^R{l)^R{l)*^R{l) 
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„2d(c)L _ 1. ^R{l)*^R{l)^L{l)*^L{l) 
e -D2 — -^'J4{A,B,X,Y)'^jAX '^lAv'^iBY '^iBX 

- 2 h{A,B,X,Y) M-- M- - Cfjx* CiAY CiBY Cfsx > (38) 

B'f'' = 0, (40) 
B^^"" = B^^%^n (a = l,---,4), (41) 



where 



■T f d Iv Jv / A n\ 

iJi[A,B,x,Y) = j (27r)4 (A;2 - )(A;2 - M|^)(A;2 - m}^){k^ ' ^IJ 

_ r d^k 1 , , 

' 4(A,B,x,y) = y (27r)4 (A;2 - M|^)(A;2 - M|^)(A;2 - m?^)(A;2 - m?^) ' ^ ^ 

Here, M^^ and m^-^ denote neutralino mass and charged slepton mass in the neutrahno 
contribution, and chargino mass and sneutrino mass in the chargino contribution. 

3.1.3 ji-e conversion in nuclei 

Finally, we give the formulas for the ix-e conversion in nuclei, ie., the process {ix+{A, Z) 
e+(A, Z)) where Z and A denote the proton and atomic numbers in a nucleus, respectively. 
The contribution again consists of the Penguin-type diagrams and the box-type diagrams. 
The box-type Feynman diagrams are depicted in Fig. 6(b) and (c). We give the effective 
Lagrangian relevant to this process at the quark level. We find that the Penguin-type 
diagrams give the following terms. 



q=u,d 

^il ^ ?l±^q^^q e^iFLPL + FnPR)^i, (44) 



g=u,a 



where the first term comes from the Penguin-type diagrams of photon exchange and the 
second one Z boson exchange. The coefficient Ql^ denotes the electric charge of the quark 
q and ^I(r) is Z boson coupling to the quark qL{R) such as Eq. (29). 
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The box-type diagrams give 

I^Mi = E er{DiPL + DfPn)i^, (45) 

q=u,d 

with 

jjL,R ^ j^{n)L,R ^ jjic)L,R ^u,d). (46) 

The coefficients are calculated to be 

-i^x°.^xO.V,B,x,y)(<i?A^e1>.'S<ir 
jj{n)R ^ D^-^%^n {q = u,d), (48) 



and 



^2r>W^ _ ^7 ^R(l)* ^R{1) ^R{d) ^R{d) 



g-4(A,S,X,y)<^^AX ^eBX^dAY^dBY 

\M^-M^-MA,B,x,y,C^T^^^ (49) 



4 Xa 

Jin^c)L _ It ^R{1)* ^R{1) ^L{u)* ^L{u) 

^ J-'u — --g'JKA,B,X,Y)'^t,Ax'^eBx'^uAY '^uBY 

+ i^xi^x.^4(AB,x,y)Ci*c2'icS^ (50) 

Note that we only take account of the vector contributions for the quark currents. 
The reason is given as follows. In the limit of the low momentum transfer which is 
appropriate for the present case {(f ~ —nT'^), we can treat the hadronic current in the 
non-relativistic limit. Furthermore, the contributions from the coherent process dominates 
over the incoherent ones if we concentrate on the relevant process such as fj, +2I Ti 
e +22 Ti. Then, the matrix element for the /i-e conversion process is dominated by the 
contribution from the vector currents. 

3.2 Decay rates and conversion rate 

Now it is straightforward to calculate the decay rates and the conversion rate, using the 
amplitudes (or the effective Lagrangian) given in the above subsection. 
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3.2.1 Ij ^ 1^ 7 decay rate 

The decay rate for Ij — > 7 is easily calculated using the amphtude (14), 



e2 



167r 

3.2.2 ir l~ If decay rate 



- C 7) - T7-<(l^2 r + (51) 



Using the expressions for the amplitude, we can calculate the decay rate. 



e 



4 



r(/,--./r /r /+) = ——ml \A[\' + \A'^\' - 2{A[Ar + A^.Af* + h.c 



where 



5127r3 



L|2 I I aR\2 r)/ /iL aR* I /iL /lii'* 



,Li2 , I .i?i2N / 16, mu 14 



+(|^^|^ + |A^|^) -1^ 



2m/^ 9 



+6(|5f |2 + |5f |2) - ^{B^B^* + 53«5f + /i.e.) 
+^{A^B^* + A^B^* + A[B^* + Af * + /i.e.) 
-^(A^5f* + A^B^* + ^^^2^* + A^B^* + h.c.) 



+1 {2(|Fii|' + \FRRf) + \Flr\^ + \Frl\^ 
+{B^Fl^ + F*^ + B^FIj, + B^F*^^ + h.c.) 
+2(Af + F*^ + h.c.) + {A[Fl^ + F*^ + /i.e.) 
-4(A^F;^ + ^^F^^ + /i.e.) - 2(A^F^^ + A^Fl^ + h.c.)] 

(52) 

(53) 
(54) 
(55) 

(56) 

Numerically, we find that a Penguin-type contribution involving A2 and dominates 
over the other contributions. In the large tan /3 region, its effect is enhanced due to the 



Fll = 


FlZ^ 




m| sin^ cos^ 


9w 


Frr = 


Fll\l^r ) 




Flr — 


FlZ'r 




m| sin^ 9w cos^ 


9w 


Frl — 


FlrIl^r- 
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same mechanism as in the case of Ij — > l~ 7 process. Furthermore, even in the case where 
tan (3 is not so large, the contribution of the Penguin-type diagram dominates over the 
box contribution, because of the logarithmic term in Eq. (52) which is quite larger than 
the other terms. ^ Then, the above formula is greatly simplified, and one finds a simple 
relation 

Brjl- ^ /-/-/+) ^ a /16 m,^. U\ 
Br{lJ ^ la) Stt 1, 3 2mi^ 9 )' ^ ' 

3.2.3 fi-e conversion rate {n + {A,Z) ^ e + {A, Z)) 

Once we know the effective Lagrangian relevant to this process at the quark level, we can 
calculate the conversion rate [7] , 

V{^,^e) = ^a^^\F{q)\'ml[\Z{A{-A^)-{2Z + N)D^^-{Z + 2N)D]i\-' 

+|Z(^f - A^) - {2Z + N)D^ -{Z + 2N)Df\^] , (58) 



where 



f)L _ j-)L : ^1 + Zr /rq\ 

^ ^^2 mlsm^ewcos^9w' ^ ' 

Df = D^Il^r iq = u,d), (60) 

and Z and N denote the proton and neutron numbers in a nucleus, respectively. Zeff has 
been determined in [6] and F(g^) is the nuclear form factor. In UTi, Z^jj = 17.6, F{q'^ ~ 
-ml) ~ 0.54 [7]. 

4 Results of the Numerical Calculations 

In this section, we present results of our numerical analysis. 

As was discussed in Section 2, we assume the universal scalar masses. Also for sim- 
plicity, we consider the so-called GUT relation among the gaugino masses 

Ml M2 M3 

— = — = — • (61) 

gf 92 93 



^This logarithmic term is obtained as a result of the phase space integration of the fermions in the 
final state, since we have an infrared singularity in the limit of nii- 0. 
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Then the SUSY breaking terms have four free parameters; the universal scalar mass (mo) , 
the SU{2)l gaugino mass at low energies (M2), the universal 74-parameter (^4 = amo) 
and mixing parameter of the two Higgs bosons {B). 

Concerning the SUSY invariant Higgs mass ^ and S-parameter which parameterize 
the mixing among hi and /i2, we determined them so that the two Higgs doublets have 
correct vacuum expectation values {hi) = f cos and (/12) = V sm P/V2. With this 

radiative electroweak symmetry breaking condition [4], we determine the mass spectra and 
mixing matrices of the supcrparticlcs. Then, wc carefully investigate the parameter space 
where tan/3 is large and masses of superparticles (especially, sleptons and electroweak 
gauginos) are quite light enough to enhance the LFV rates. As a result, we found that 
there indeed exists parameter space where the above conditions are satisfied. We checked, 
for M2 — 80GeV, tan /3 can be as large as about 50.^ This result implies that there are 
regions in the parameter space where the LFV processes have large branching ratios due 
to the large tan /3 enhancement mechanism.^ 

We also put constraints from experiments. Besides our requirement that the light- 
est superparticle be neutral, we use consequences of the negative searches for the su- 
perparticles [8]. Wc also impose a constraint on SUSY contribution to the anomalous 
magnetic dipole- moment of the muon [11, 12]. The experimental value of ^{g — 2) is 
1165923(8.4) x 10^^ [8]. On the other hand, the theoretical prediction of standard model 
is 11659180(15.3) x 10"^° or 11659183(7.6) x 10"^° [12], where the difference is due to 
different estimates of hadronic contributions. In our paper, we adopt the first one in order 
to derive conservative bound. Therefore, the SUSY contribution should be constrained as 

- 26.7 X 10"^ <{g- 2)f^^ < 46.7 x 10"^ (62) 

where two sigma experimental error is considered. The SUSY contribution is shown in 
Fig. 7. Here, we take the parameter a = at the gravitational scale and M2 = 100 GeV 
at low energies. The horizontal line is taken to be the left-handed selectron mass with 

^Throughout this paper, we take the top quark mass nit = 174 GeV [8]. Also we take the bottom 
quark mass nib = 4.25GeV [9], which corresponds to 3.1 GeV at the Z mass scale. 

^Note that the situation here contrasts to the case of the Yukawa unification where the radiative 
breaking with the universal scalar mass requires heavy superparticle spectrum, larger than, say, 500 GeV 
[10]. 
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the D-term contribution, which we denote by mg^. One finds that a significant region 
of the parameter space is excluded by this constraint in the large tan/3 region. This is 
because the same enhancement mechanism as the LFV processes works in the diagrams 
contributing to the g — 2. For completeness, we will give formulas of the contribution of 
the superparticle loops to the anomalous magnetic dipole- moment in Appendix C. 

Let us now discuss the branching ratios for each LFV process. First we consider the 
case where the neutrino mixing matrix is described by the KM matrix. 



As the first trial, we shall consider the case where V — Vkm, where we take S12 = 0.22, 
S23 = 0.04 and S13 = 0.0035 in the standard notation [8]. We ignore the possible 
Kobayashi-Maskawa complex phase and consider V to be real, for simplicity. The eigen- 
values of the neutrino Yukawa couplings are assumed to be equal to those of the up-type 
quarks at the gravitational scale. Since the magnitude of the top quark Yukawa coupling 
is close to its perturbative bound, this ansatz will maximize the magnitude of LFV in the 
slepton mass matrix. Also, to determine the right-handed neutrino Majorana mass Mr, 
we fix the tau neutrino mass at 10 eV so that it constitutes the hot component of the 
dark matter of the Universe. In this case, Mr is about 10^^~^^ GeV. 

Solving the RGEs numerically, we obtain the mass squared matrix for the SU{2)l 
doublet sleptons at the electroweak scale 



where tan (3 varies from 3 to 30, M2 = and a = 0. For a non- vanishing M2, the diagonal 
elements of the above matrix become larger and the fiavor- violating off-diagonal elements 
become relatively less important, as the gaugino mass gets larger. Effect of non- vanishing 
a-parameter can be seen from Eq. (7), which does not change the result drastically. In 
the following numerical calculations we will take a = 0. 

We find in Eq. (63) the off-diagonal elements in the mass matrix are small. This is 
because the off-diagonal slepton masses are proportional to V^^V^^ in the case of hierar- 



4.1 Case 1) V = Vkm 




f 1.00 

(0.30 - 0.43) X 10-^ 
V -(0.74 - 1.07) X 10-^ 



(0.30 - 0.43) X 10-4 
1.00 

-(0.54 - 0.78) X 10-2 



(0.74 - 1.07) X 10-^ \ 
(0.54 - 0.78) X 10-2 
0.77 - 0.80 / 




(63) 
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chical neutrino masses, which are small if we assume that V is equal to the KM matrix. 
Nevertheless, as will be shown shortly, the enhancement in the large tan /3 region yields 
large branching ratios for the LFV processes, which is close to the present experimental 
upper bounds. 

4.1.1 — > 67 

Result of our computation on the branching ratio Br{fi 67) is shown in Fig. 8 for 
M2 = 100 GeV. The horizontal line is taken to be the left-handed selectron mass with 
the D-term contribution, mg^. Real lines are for > 0, while dashed lines for /j, < 0. 
We can find that the branching ratios are rather insensitive to the choice of the sign of 
the /x-parameter, in particular when tan (3 is large. For the large tan (3 case, some regions 
of small slepton masses are excluded by the constraint from g — 2. As can be seen from 
Fig. 7, it is less stringent for // > case than jj, < caseJ One can see that even if we 
impose this constraint, the branching ratio can be as large as 10^^^, which is very close 
to the present experimental bound -Br(/i — > e7)|exp < 4.9 x 10"^^. For smaller value of 
tan /5, the branching ratio reduces obeying oc tan^ /3. 

We compared the chargino loop contribution with the neutralino loop contribution 
and found that the former dominates. This is important when we compare our results 
with the case of SU{5) grand unification. (See Section 5.) 

In Fig. 9, we show the case of M2 — 200 GeV. The maximum of the branching ratio 
is about 10~^^ for tan/? = 30, about one order of magnitude smaller than the M2 — 100 
GeV case. We also studied the case M2 — 80 GeV, and found that the branching ratio is 
about factor 2 larger than the M2 — 100 GeV case. 

''Here, we should comment that the SUSY contribution to the 6 — > 57 process is also significant and 

some part of the parameter space should be excluded [12, 13, 14]. However, it is complicated to estimate 
the SUSY contribution to the b ^ process, since the chargino loop can contribute either constructively 
or destructively to the others, especially charged Higgs boson loop. Thus, it seems to us that to determine 
which regions of the parameter space are really eliminated contains some delicate issues as discussed by 
Ref. [14]. We believe that such an analysis is out of the scope of our paper, but a work in a future 
communication. Thus, we do not use the constraint from the 6 ^ 57 process. 
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Next, let us consider the process fi~ e~e~e^ . Currently the experimental upper bound 
on the branching ratio of this process is 1.0 x 10^^^ [8]. We show results of our calculation 
to this process in Fig. 10 for M2 = 100 GeV. The branching ratio has the maximum of 
~ 10^^'^ for the large tan/? with the small gaugino mass. One can check that this process 
is dominated by the Penguin-type diagrams. Indeed compared with the branching ratio 
of — > 67, one finds a simple relation 

14^^^ - r X 10-3, (64) 

which is in agreement with the ratio expected by the dominance of the Penguin-type 
diagrams, Eq. (57). 

4.1.3 iJ,-e conversion in ffTi 

Experimentally, //-e conversion rate in nuclei is also constrained strongly. The experi- 
mental upper bound on the conversion rate with the target g^Ti reaches 4.3x10"^^ [8]. 
We show results of our calculation to this process in Fig. 11 for M2 — 100 GeV. The 
branching ratio takes its maximal value of ~ 10"^^ in the parameter region where tan /3 
is large and the gaugino masses are small. On the other hand, for the small tan/3 and 
< the cancelation among the diagrams occurs and the event rate damps rapidly. The 
Penguin-type diagram is not dominant in the small tan (3 region because there is not the 
same logarithmic enhancement as /x~ — > e~e~e'^. 

4.1.4 T ^ //7 

Finally we would present our result for r — > /X7 in Fig. 12. We find with M2 — 100 GeV, 
the branching ratio is as large as 10~^, one and a half order of magnitude smaller than 
the present experimental bound Sr(r — > /X7)|exp < 4.2 x 10~^ [8]. Similar to the case of 
II 67, it can be seen that the branching ratio is proportional to tan /3 squared. 
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4.2 Case 2) Neutrino mixing implied by atmospheric neutrino 
deficit 



A class of solutions to the atmospheric and solar neutrino deficits requires a maximal 
mixing of the tau and muon neutrinos, yielding a large off-diagonal element in the slepton 
mass matrix. The neutrino mixing matrix we take in this example is 



and the tau neutrino mass is assumed to be 0.4 eV [15]. Here, we only consider the 
generation mixing of the second and third generations and ignore the others. The (1,3) and 
(3,1) elements of the mixing matrix cannot be determined from the solar and atmospheric 
neutrino deficits. This uncertainty, however, does not matter if we only consider the 
LFV process among the second and third generations. As in the case 1), wc assume the 
magnitude of the third generation neutrino Yukawa coupling f^s is equal to the top quark 
Yukawa coupling at the gravitational scale. The latter choice will give us a maximum 
violation of LFV in the slepton mass matrix. 

The result for Sr(r — > /X7) is shown in Fig. 13. We find that in some portion of the 
parameter space, the branching ratio exceeds the present experimental upper bound, in 
particular when tan /3 is large and the superparticles are light. 

5 Conclusions and Discussion 

In this paper, we have considered LFV in the minimal supersymmetric standard model 
(MSSM) with the right-handed neutrino multiplets. In the presence of the Yukawa cou- 
plings of the right-handed neutrinos, the left-handed slepton mass matrix, m|, loses its 
universal property even if we assume the minimal supcrgravity type boundary condition 
on sfermion masses. In our case, due to the renormalization effect, as can be seen from 
Eq. (63), wc obtain LFV in m| as well as smaller value of (3,3) clement of m| compared 
with the other diagonal elements, which is typical feature of the case with right-handed 
neutrino [16]. We have calculated the interaction rates for the various LFV processes with 
the full diagonalization of the slepton mass matrices and of the chargino and neutralino 




0.87 X 10" 
0.755 
-0.656 



0.656 
0.755 



(65) 
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mass matrices. We emphasized the enhancement of the interaction rates for large tan (3, 
the ratio of the VEVs of the two Higgs doublets. This enhancement is originated to the 
fact that there is a freedom to pick up one of two vacuum expectation values in the MSSM 
in the magnetic dipole-moment type diagrams. For example, for the process Ij — > li'j, the 
diagrams of the type Fig. 1(c) and Fig. 2 give the enhancement. Even when the mix- 
ing matrix in the lepton sector has a similar structure as the KM-matrix of the quark 
sector, the enhancement mechanism can make the branching ratios close to the present 
experimental bounds. 

It is interesting to compare the LFV processes induced by the right-handed neutrino 
Yukawa couplings with those in the minimal SU{5) grand unified theory [17, 18]. In the 
latter case, the renormalization-group flow above the GUT scale results in LFV in the 
SU (2) L singlet (right-handed) slepton masses. Let us consider, for example, the resulting 
branching ratio of — > 67. The diagrams which will give the enhancement in the large 
tan/3 region are similar to Fig. 1(c) and 2(a). The important difference from the previous 
case is on Fig. 2. Now, only the diagrams involving the bino contributes, since the wino 
does not couple to the singlet sleptons. In this case, we can see that contributions coming 
from the two diagrams Fig. 1(c) and Fig. 2(a) have opposite signs, and thus partially 
cancel out with each other. Numerical result is shown in Fig. 14. The horizontal line is 
the mass of the right-handed selectron with the D-term contribution, mg^. Here, we have 
taken M2 — 100 GeV. The branching ratio never exceeds 10~^^, more than two orders of 
magnitude beneath the present experimental upper bound. Also one flnds regions where 
the branching ratio becomes very small due to the cancelation explained above. 

What happens if the standard model with the right-handed neutrinos is embedded 
in the framework of SU (5) GUT? In this case, both the mass matrix of the left-handed 
sleptons and that of the right-handed ones have LFV. The situation is quite similar to 
the case of 150(10) GUT [18, 19]. For example, if we consider the /i ej, the dominant 
diagram will be similar to Fig. 1(c), which however picks up {'m\^)f^, proportional to tau- 
lepton mass. Thus we expect further enhancement in the branching ratio by {m-j-/mfj)'^ 
compared to the case we studied in this paper. 

To conclude our paper, we should emphasize that the branching ratios of the LFV 
processes induced by the right-handed neutrino Yukawa couplings can be close to the 
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present experimental bounds and can be within the reach of future experiments. Efforts 
of searching for these LFV signals should be encouraged. 
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A Renormalization Group Equations 



In this appendix, we give the one-loop renormalization group equations (RGEs) for the 
Yukawa couplings and the soft SUSY breaking terms in the scalar potential. The RGEs for 
the gauge coupling constants and the gaugino masses are unchanged at the one-loop level, 
since the right-handed neutrinos are singlet under the standard-model gauge symmetry. 

• Yukawa coupling constants 

1 rr 9 

T67r2 . 



--gl - 'igl + 3Tr(/,/i) + Tr(/j/) \ //^ 



167r2 



:gt - 2>gi + 2,Tr{Ul) + Tr(/,/J) fj 



t f\i3 



(66) 



(67) 



Soft breaking terms 



167r2 

+2 {flmlfi + rhljIfi + AlAil 



6 



- ( ^-gl \M,f + 6^i iM^l'j 5i - -glS5{ 



167r2 



2 (m?/j/ + /,//m?); 



(68) 
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where 



+4 (/,m|//+mU// + M)'. 



2{mlfJl + fJlmiy 



167r2 

+4 (Umlfl + ml^Ufl + A^Ai) 
1 



-pi - ^gl + 3Tr(/]/,) + Tr(///,) ^ 4 



+2 [--glMr - 3glM2 + 3Tr(/J^,) + Tr(//A0} Z/^' 



-pi - ^gl + 3rr(/t/„) + TrUlU) \ A^i 



167r2 



(69) 



(70) 



+2 {--^i^Mi - 3^2M2 + 3rr(/t^„) + rr(/t^.)} fj 
+4(/,/tA,)^^- + 5(A,/t/,)^^- + 2{U}A,r + (^.///O^'l , (71) 



S = Tr{m\ + m| — 2mu — m\ + m?) — m^^ + m^2- 



(72) 



Here, we followed the GUT convention for the normalization of U{1)y gauge coupling 
constant gi, such as gy = ^gf- 



B Interaction of gaugino-sfermion-fermion 

In this appendix, we give our notations and conventions adopted in Section 3 and give 
vertices relevant for our calculation. 

Let us first discuss fermions. We denote by li, Ui and di the fermion mass eigenstates 
with the obvious meaning. The subscript i (i = 1, 2, 3) represents the generation. As for 
the neutrinos, their masses are small and negligible. In our convention, is the SU{2)l 
iso doublet partner to eu- 

Next we consider sfermions. Let fu and fm be the superpartners of fu and fm, 
respectively. Here, / stands for Z, u or d. The mass matrix for the sfermions can be 
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written in the following form, 

{firn)(^^2^^ If ) (£)' (^3) 

where m| and m|j are 3x3 hermitian matrices and m|,^ is a 3 x 3 matrix. These elements 
are given from Eqs. (1,3) as following, 

ml = m)^ +mj + m| cos 2/3(7;^^ -gf^sin^M, (74) 

^ "^/^ + "4-"^|cos2/3(T3^^-g{,„sin2^M/), (75) 

2 ^ f - Af V sin p / - nifH cot f3 (f = u), 

\ AfV cos P/V^-ruf iitanP {f ^ d,l), ^ ^ 

where T^j^^p,-^ and are weak isospin and electric charge respectively. Here, = 



for squarks, = m? for sleptons, and m ^ arc each right-handed sfermion soft-breaking 
masses. We assume the above mass matrix to be real. This is, in general, not diagonal 
and include mixing between different generations. We diagonalize the mass matrix M."^ 
by a 6 X 6 real orthogonal matrix as 

M^U^^ = (diagonal) , (77) 

and we denote its eigenvalues by (X = 1, • • • , 6). The mass eigenstate is then written 
as 

fx = Uy~u + ^i,+3^, (X = 1, ■ ■ ■ , 6). (78) 

Conversely, we have 

fu= Ugfx =UUx: (79) 
fm= U,Z,xfx =U^,+sfx. (80) 

An attention should be paid to the neutrinos since there is no right-handed sneutrino 
in the MSSM. Let uu be the supcrpartner of the neutrino Ui. The mass eigenstate ux 

{X = 1, 2, 3) is related to uu as 

i>u = U^xi^x. (81) 
We now turn to charginos. The mass matrix of the charginos is given by 

M2 V2mwC0sP \ f Wj^ 



- ( H,,) ^ ~ + .... (82) 
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This matrix Mc is diagonalized by 2 x 2 real orthogonal matrices Ol and Or as 

OrMcOI = (diagonal). (83) 



Define 



Then 



XlL 
X2L 



Oj 



H 



IL 



Xa — Xal + Xar 



forms a Dirac fermion with mass M--. 



X2R J 

(A =1,2) 



(84) 
(85) 



Finally we consider neutralinos. The mass matrix of the neutralino sector is given by 

( Bl \ 



h.c, 



(86) 



where 



M, 



N 



( M\ — sin cos/? sin sin /9 \ 

M2 cos cos /? —mz cos sin /? 

—raz sin 6*14/ cos (5 mz cos Qw cos /3 — // 

V sin sin /? — cos Qw sin /3 — 



The diagonalization is done by a real orthogonal matrix O^, 

OjsiMjsiO\; — diagonal. 



The mass eigenstates are given by 
where 



(A,S = l,---,4) 



= (-Sl, W^l, HIl)- 

We have thus Majorana spinors 

x\ = x\l^x\r. (A = 1,---,4) 



/ 

(87) 



(88) 



(89) 



(90) 



(91) 



with mass M^^o . 
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We now give the interaction Lagrangian of fermion-sfermion-chargino, 

Ant = k{C^AxPR + C^AxPl)Xa'^X 

+d,icS'JPR + Cf'}'M)XAUx 

+u,{cSSPn + C^SPL)x\dx + h.c, (92) 

where the coefficients are 

CTx = 92 r- J Ql)a2^^„ 

V 2mvK cos p 

= 5'2^^^ ^(OL)A2C/j(:,i+3, 

V2mvi/ cos p 

= g2{-(Qi^)Al^l.+ J Qil)A2^W^ 

V 2mty sm p 

^iAX - 92—?^ -yJL)A2'Jx,v 

y IVTiy/ COS p 

GfS = g2{-(QL)Al^i,+ ^'"'^^ J QL)A2^i,+3}, 

v2mvKCOsp 

Cltl = 92 A 0n)A2Ui,. (93) 

V 2mvt^ sm p 

The interaction Lagrangian of fermion-sfermion-neutrahno is similarly written as 

Ant = UN^i^^PR + NfigPL)xVfx (94) 
where / stands for /, i/, c? and u. The coefficients are 



fl'i = -^{[-{ON)A2-{0^)Aitanew]U^x, + ^^iO^)AsU^x 
V 2 cos p 



= -^{-^^(0^)^3C/^, + 2(0^)^itane;^C/^,^3}, 

V z TTivy cos p 

= -;^[(0iv)A2-(0^^)Altan^^]t/^,„ 

= -;^{[-(0iv)^2 + ^(0iv)Aitan^H.][/i,, + -^^( 
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V 2 niw cos /? '3 

N-S = -^{^^A0N)A.Ul.-\i^^9w{0M)AiUl,^,}. (95) 

V 2 sm p '3 ' 

C Anomalous magnetic dipole-moment of the muon 

The magnetic dipole-moment interaction of muon is written as the following form; 

F{q')u{pf)a,,q^t''u{pi), (96) 



2m^ 



where q — Pf—Pi and e the polarization vector of external photon. Then, the anomalous 
magnetic dipole-moment of muon is 

(^-2)^ = 2F(?2 = 0). (97) 

We can write SUSY contributions as {g - 2)^^-^^ = {g^'^'i + g'^^^)^i- The first term gj^^ 
represents the chargino-loop contribution as 

(C) _ 1 ir^m 2 + 3a; AX - Gx\x + ^^ax + ^^AX In a: ax 

^'^ - 487r2 mij^^^^^l (1 - xax)^ 

_^ 1 '^i^^xa r<L(i) ^R(i)* -3 + 4xAx - x\x - 2 In xax 

+ 167r^ mi, ^2Ax^2Ax (1 - ^^^)3 

+(L ^ i?) (98) 

where x^x = / . 

The neutrahno-loop contribution g'^^^ is 

(jv) 1 TTT'l I ^r^(0 |2 1 - 6a; AX + 3a:^^ + 2a:^;^ - 6x\x In xax 

~ 487r2m? (1 - xax)^ 

1 ^^L(i) ^rR(iy 1 - a^ix + 2.TAX In xax 

IGtt^ m?^ ^"^^ '^^'^ (1 — xax)'^ 

+ (L ^ i?), (99) 
where a; ax = /m3 . 

Xa' i-x 
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Figure 1: Fcynman diagrams which give rise to Ij li^j. The symbols en, uu- B, W^, and 
W~ represent left-handed charged sleptons, left-handed sneutrinos, bino, neutral wino, 
and charged wino, respectively. In (a) and (b), the blob in the the slepton/sneutrino line 
indicates the flavor-violating mass insertion of the left-handed slepton and another blob 
in the external line the chirality flip of the external lepton Ij. In (c), the blobs in the 
slepton line indicate the insertions of the flavor- violating mass ijnYi) ^"^^ left-right 
mixing mass (^iijjj), and another blob in the bino line the chirality flip of the bino B. 

Figure 2: Feynman diagrams which give rise to the large tan/3 enhancement due to the 
gaugino-higgsino mixing in the process of Ij li^y. The symbols en, Un, B, W^, W~, 
H^, and H represent left-handed charged sleptons, left-handed sneutrinos, bino, neutral 
wino, charged wino, neutral higgsino, and charged higgsino, respectively. The blob in the 
slepton/sneutrino line indicates the insertion of the flavor- violating mass (mil). The blobs 
in the gaugino-higgsino line indicate the mass insertions for gaugino-higgsino mixing, that 
is, /i denotes higgsino {H1-H2) mass mixing, v sin j3 the gaugino-higgsino {H2-W) mass 
mixing, and M2 the wino mass. The value of tan (3 comes from Yukawa coupling constant 
fi. ~ 1/ cos (3 and v.e.v. of /i2, f sin/5. 

Figure 3: Feynman diagrams for the process Ij li'-f. (a) represents the contributions 
from neutralino Xa and slepton Ix loop, and (b) the contributions from chargino Xa 
sneutrino ux loop. 

Figure 4: Penguin-type diagrams for the process /~ KK^t which photon 7 and 
Z-boson are exchanged. The blob indicates Ij-k-j vertex such as Fig. 3 or Ij-li-Z vertex 
where Z-boson is external. 

Figure 5: Box-type diagrams for the process Ij l^l^lf . Here, (a) represents the 
contributions from neutralino Xa and slepton Ix loop, while (b) the contributions from 
chargino Xa sneutrino vx loop. 

Figure 6: Feynman diagrams for the process ji-e conversion at the quark level. In (a), 
the Penguin-type diagram is depicted. The blob indicates Ij-k-j vertex such as Fig. 3 
or Ij-li-Z vertex such as Fig. 4. In (b) and (c), the box- type diagrams are depicted; i.e., 
(b) represents the contributions from neutralino x% slepton Ix and squark qx {q = u, d) 
loop, and (c) the contributions from chargino x~^, sneutrino vx and squark qx {q — u, d) 
loop. 
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Figure 7: The values of the SUSY contribution to the anomalous magnetic dipole-moment 
of muon {g — 2)^^^^ as a function of the left-handed selectron mass with the D-term 
contribution, which we denote by mg^. Here we assume a = at the gravitational scale. 
Real lines correspond to the case for fi > 0, while dashed lines for /i < 0. Here we have 
taken M2 = 100 GeV and tan/3 = 3, 10, 30. The shaded regions are excluded by the 
present experiments. 

Figure 8: Branching ratios for the process /x — > 67 in the Case 1) V = Vkm as a function of 
the left-handed selectron mass with the D-term contribution, mg^ . Real lines correspond 
to the case for > 0, while dashed lines for < 0. Here we have taken M2 = 100 GeV 
and tan/3 = 3, 10, 30. We also show the present experimental upper bound for this 
process by the dash-dotted line. 

Figure 9: Same as Fig. 8 except for M2 = 200 GeV. 

Figure 10: Branching ratios for the process /i~ e~e~e^ in the Case 1) V = Vkm as 
a function of the left-handed selectron mass with the D-term contribution, mg^. Real 
lines correspond to the case for > 0, while dashed hues for // < 0. Here we have taken 
M2 = 100 GeV and tan (3 — ?>, 10, 30. We also show the present experimental upper 
bound for this process by the dash-dotted hue. 

Figure 11: The /x-e conversion rates in nuclei ^fTi in the Case 1)V = Vkm as a function of 
the left-handed selectron mass with the D-term contribution, mg^. Real lines correspond 
to the case for /i > 0, while dashed lines for fi < 0. Here we have taken M2 = 100 GeV 
and tan/3 = 3, 10, 30. We also show the present experimental upper bound for this 
process by the dash-dotted hue. 

Figure 12: Branching ratios for the process r — > /i7 in the Case 1) V — Vkm as a 
function of the left-handed selectron mass with the D-term contribution, mg^. Real lines 
correspond to the case for /i > 0, while dashed lines for /x < 0. Here we have taken 
M2 = 100 GeV and tan/3 = 3, 10, 30. We also show the present experimental upper 
bound for this process by the dash-dotted line. 

Figure 13: Branching ratios for the process r — > /X7 in the Case 2) neutrino mixing 
implied by atmospheric neutrino deficit, as a function of the left-handed selectron mass 
with the D-term contribution, mg^. Real lines correspond to the case for > 0, while 
dashed lines for < 0. Here we have taken M2 = 100 GeV and tan/3 = 3, 10, 30. We 
also show the present experimental upper bound for this process by the dash-dotted line. 
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Figure 14: Branching ratios for the process // — > 67 in the case for the minimal SU{5) 
grand unified theory, as a function of the right-handed selectron mass with the D-term 
contribution, mg^. Here we have taken // > 0, M2 = 100 GeV, and tan/? = 3, 10, 30. We 
also show the present experimental upper bound for this process by the dash-dotted line. 
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